Statistics of resonances and delay times: A criterion for Metal-Insulator transitions 



(N 
O 

o 



Tsampikos Kottos 1 , and Matthias Weiss 2 
1 Max-Planck-Institut fur Stromungsforschung, Bunsenstrafie 10, D-37073 Gottingen, Germany 
2 EMBL, Meyerhofstr. 1, D-69117 Heidelberg, Germany 

We study the distributions of the normalized resonance widths V(T) and delay times Vij) for 3D 
disordered tight-binding systems at the metal-insulator transition (MIT) by attaching leads to the 
boundary sites. Both distributions are scale invariant, independent of the microscopic details of the 
random potential, and the number of channels. Theoretical considerations suggest the existence of a 
scaling theory for "P(r) in finite samples, and numerical calculations confirm this hypothesis. Based 
on this, we give a new criterion for the determination and analysis of the MIT. 
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Quantum mechanical scattering has been a subject of 
intensive research activity during the last years 
both in Mesoscopic Physics and in studies of Quantum 
Chaos. Among the most interesting quantities for the 
description of a scattering process are the Wigner de- 
lay times and resonance widths. The former quantity 
captures the time- dependent aspects of quantum scat- 
tering. It can be interpreted as the typical time an al- 
most monochromatic wave packet remains in the inter- 
action region. It is related to the energy derivative of the 
total phase shift $(£') = — ilndet S(E) of the scatter- 
ing matrix S(E), i.e. t(E) = jj d ^ d ^ , where M is the 
number of channels. The resonances represent interme- 
diate states with finite lifetimes to which bound states 
of a closed system are converted due to the coupling to 
a continuum. On a formal level, resonances show up 
as poles of the S'-matrix occurring at complex energies 



f — E 



; r„, where E n and Y n are called position and 



width of the resonance, respectively. 

For chaotic/ballistic systems many results are known 
concerning the distributions of resonance widths V(T) 
and Wigner delay times V(t) Recently, the interest 
has extended to systems showing diffusion and localiza- 
tion, where V(T) j|§ and V(t) j| were found to follow 
universal distributions with algebraic decay. At the same 
time, an attempt to understand systems at critical con- 
ditions, was done in Q where V{T) and V(t) for one- 
dimensional (ID) quasi-periodic systems was studied. 

Despite the progress in understanding the scattering 
process of various systems, a significant class was left out 
of the investigation. These are systems, whose closed 
analogues show a metal-insulator transition (MIT) as an 
external parameter changes. In the metallic regime, the 
eigenstates are extended and the statistical properties of 
the spectrum are quite well described by random matrix 
theory In particular, the level spacing distribution 
is very well fitted by the Wigner surmise. Deep in the 
localized regime, the levels become uncorrelated leading 
to a Poissonian level spacing distribution, and the eigen- 
functions are exponentially localized. At the MIT the 
eigenfunctions exhibit multifractal behavior with strong 
fluctuations on all scales while the eigenvalue statistics is 



characterized by a third universal distribution |7]||. At 
the same time a considerable effort was made to under- 
stand the shape of the conductance distribution V(g) at 
MIT d |ll[. However, it is still unclear whether the lim- 
iting V(g) is entirely universal, i.e. dependent only on 
the dimensionality and symmetry class, as required by 
the one-parameter scaling theory of localization Q . The 
latter is one of the major achievements in the long his- 
tory of studying the MIT. Its basic assumption is that 
close to the MIT the change of the conductance g with 
the sample size L depends only on the conductance itself, 
and not separately on energy, disorder, size and shape of 
the sample, the mean free path etc. 

The most prominent realization of systems which un- 
dergo a MIT (with increasing strength of disorder) is the 
three-dimensional (3D) Anderson model ||. Although 
a lot of studies have been devoted to the analysis of 
eigenfunctions and eigenvalues @-|| and of conductance 
]l0| , pl| at the MIT, the properties of delay times and res- 
onances of the S'-matrix were left unexplored. 

Here, for the first time we address this issue and 
present consequences of the MIT on the statistical prop- 
erties of the (properly) rescaled resonance widths T and 
delay times f. We find that they follow a new univer- 
sal distribution, i.e. independent of the microscopic de- 
tails of the random potential, and number of channels 
M jiJI . Specifically, they decay asymptotically with pow- 
ers which are different from those found for diffusive or 
localized systems i.e. 
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r = r/A 

f = tMA 



(i) 



where A is the mean level spacing. We relate the power- 
law decay of V(T) to the anomalous diffusion at the MIT. 
Finally, based on the resonance width distribution, we 
suggest a new method for determining and analyzing the 
emergence of the MIT and propose a scaling theory near 
the critical point. Our results are confirmed by numerical 
calculations and are supported by theoretical arguments. 

We consider a 3D sample of volume L 3 . To each site of 
the layer n x = 1 we attach M = L 2 semi-infinite single- 
mode leads, the simplest possible multichannel scattering 
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set up |l]|. The mathematical model that describes the 
sample is the tight-binding Hamiltonian (TBH) 



< T > = 



H L =Y,\n)V n (n\+ £ |n)(m| 



(2) 



(n,m) 



where n = (n x ,n y ,n z ) labels all the N = L 3 sites of 
a cubic lattice, while the second sum is taken over all 
nearest-neighbor pairs (n, m) on the lattice. The on-site 
potential V n for 1 < n x , n y , n z < L is independently and 
identically distributed with probability V(V n ). We study 
three different distributions for the random potential: (a) 
A box distribution, i.e. the V n are uniformly distributed 
on the interval [— V/2, V/2], (b) a Gaussian distribution 
with zero mean and variance V 2 /12, and (c) a Cauchy 
distribution V(V n ) — ^/ V Y +V i) ■ For the system defined 
by Eq. (§) the MIT for e"~ occurs for V = V C with (a) 
V c ~ 16.5, (b) V c ~ 21.3, and (c) V c ~ 4.26 @. Each 
lead is described by a ID semi-infinite TBH 

— oo 

h m = Y J (\ n><n+1 \ + \ n+1><n \) ■ ( 3 ) 

n=l 

Using standard methods one can write the scattering 
matrix in the form Il4f 



S(E) =l-2i sin(fc) W T (E - H cS ) W 



(4) 



where 1 is the M x M unit matrix, k — arccos(i^/2) is 
the wave vector supported in the leads while 7i c ff is an 
effective non-hcrmitian Hamiltonian given by 



'HcS = Hl 



e lk WW T . 



(5) 



Here W is a N x M matrix with elements equal to zero 
or unity. It describes at which site of the sample we 
attach the leads. Moreover, since arccos(£'/2) changes 
only slightly in the center of the band, we put E = and 
neglect the energy dependence of H B s- The poles of S 
are then equal to the complex eigenvalues £ of TL c s. 
The Wigner delay time r = xf d ^ , can be written 

as 



r(E) = ^TrQ(E) ; Q(E) = - iS ^E)^-, (6) 

where Q(E) is called Wigner-Smith matrix. Using 
Eq. (0), we calculated the M x M matrix Q(E) 

Q(E) = (1 - i sin(fc) K)' 1 (cot(fc)K + 2 sin(fc) W^GEQW) 
(l + ^sin(fc)i^)- 1 , (7) 

where K = W^GW, Q = (EE - Hl)- 1 and £ = 
1 + 7}WW^ . Relation (j^) is very convenient for numerical 
calculations since it anticipates the numerical differenti- 
ation which is a rather unstable operation. The trace of 
Q is associated with the density of states of the closed 
system M and the mean Wigner delay is given by 
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(8) 



In order to investigate the distributions of resonance 
widths r and delay times f , we consider the integrated 
distribution 



Pint (a) 



V(x')dx' 



(9) 



whose derivative V(x) = —dV int /dx determines the prob- 
ability density of resonance widths V(x — T) and delay 
times V(x = f), respectively. The maximum size of the 
matrices that we used in our analysis was N = 8000 
(L = 20). For better statistics a considerable number of 
different disorder realizations has been used. In all cases 
we had at least 10000 data for statistical processing. 
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L=10 (box) 
L=14 (box) 
L=18 (box) 
L=20 (box) 
L=10 (Gauss) 
L=14 (Gauss) 
L=10 (Cauchy) 
L=14 (Cauchy) 
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FIG. 1. (a) Universal behavior of V(t) at the MIT (re- 
ported here as VQn(T))) for various sample sizes L and po- 
tential distributions, (b) The integrated distribution Vint(T) 
asymptotically decays as 'Pint(r) ~ f -0-333 ^agjjgj i me ). 

Figure |l^i displays our results for the distribution of the 
logarithm of the rescaled resonance widths V(ln(T )) |l[| 
for the three different distributions V(V n ) of the random 
potential and for various sample sizes L. The body of the 
distribution function in all cases coincides and does not 
change its shape or width. Of course, the far tail of this 
universal distribution develops better with increasing L. 
The sharp peak appearing at the right is an artifact of our 
choice to neglect the energy dependence of H e s- Taking 
the energy dependence into account results in a much 
milder behavior (see for example Figs. 4,6 in rcf. We 
conclude therefore that at the MIT, the distribution of 
rescaled resonances is indeed scale- invariant independent 
of the microscopic details of the potential. 

Next we turn to the quantitative analysis of P(T)- To 
this end we analyze numerically the integrated distribu- 
tion 'Pi n t(f ). Our results are reported in Fig. An in- 
verse power law Vint(T) ~ T~ Q is evident. The best fit to 
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the numerical data yields a = 0.333±0.005 in accordance 
with Eq. The behavior of the statistically insignifi- 
cant extreme large T tails of V(T) (associated with the 
sharp peak on the right) is essentially determined by the 
coupling to the leads, and therefore is model-dependent. 
Nevertheless, it is reasonable to assume that the relative 
number of resonances contributing to the peak is pro- 
portional to M/L 3 ~ L^ 1 and thus these extreme non- 
universal far tails subside as L increases in agreement 
with our data (see Fig. [l]) 




FIG. 2. (a) Vi nt (V, L) as a function of V for different system 
sizes L provides a means to determine the critical point V c of 
the MIT (vertical line at V = 16.5). (b) The one-parameter 
scaling of V^iV, L) [Eq. (|l3|)] is confirmed for various system 
sizes L and disorder strengths V using the box distribution. 



The following heuristic argument provides some un- 
derstanding of the universal character of V(T) and gives 
a quantitative estimate for the power law (fy). At this 
point we wish to generalize our analysis to d-dimensional 
systems with a MIT. The specific case of d = 3 will be 
recovered at the end. At V = V c the conductance of a d- 
dimensional sample has a finite value g c ~ 1. Approach- 
ing the MIT from the metallic side one has g ~ Et/A, 
where Et = D /R 2 is the Thouless energy, D is the diffu- 
sion coefficient, and A ~ l/R d is the mean level spacing 
in a d-dimensional sample with linear size R. This yields 
D ~ g c /R d ~ 2 at V c . Taking into account that D = R 2 /t, 
we get for the spreading of an excitation at the MIT 



R d (t) ~ g c t 



(10) 



We now consider the influence of the leads on the spec- 
trum of an isolated sample (for simplicity we consider 
that the leads are attached to the sites of the layer n x = 1, 
but the argument can be generalized easily). Due to the 
leads the levels are broadened by a width T and the eigen- 
states are turned into resonant states. The inverse of 
r represents the quantum lifetime of an electron in the 
corresponding resonant state before it escapes out into 
the leads. We suppose that the resonant states are uni- 
formly distributed along the sample. Their lifetime can 
be associated with the typical time ~ l/r# ~ R d / 'g c 
[Eq. din)] a particle needs to reach the boundary n x = 1, 



when starting at distance R apart. Moreover we assume 
that the particles are uniformly distributed inside the 
sample. This classical picture can be justified for all 
states with Tr > Tl — 1/£l = L d /g c ~ A since the 
level discreteness is unimportant. Therefore the number 
of states with r > Tr is 

Pint(r«) w RL d - x lL d = R/L ~ (g c /T R ) 1 / d /L . (11) 

By repeating the same argument for R — > L (i.e. Fr — » 
Tl) wc finally obtain for all T > Tl 

^i„t(r) ~ g 1 J d {\/T)V d V(f ) ~ gUdf-d+i/d) (12) 



which is equivalent to Eq. 



for the case d = 3. Since 
[ |l0| , we expect V(T) to be 



g c is universal at the MIT 
universal as well. An additional conclusion is that the 
power-law decay of V(T) is determined by the anomalous 
diffusion at the MIT, which leads to the dependence on 
the dimensionality of the system. 

In the original proposal of the scaling theory of lo- 
calization, the conductance g is the relevant parameter 
H . A manifestation of this statement is seen in Eq. jl^ ) 
where V® n t = Pint (Jo) is proportional to the conductance 
g. It is therefore natural to expect that V® nt will follow a 
scaling behavior for finite L (and for some To ~ 1), that 
is similar to the one obeyed by the conductance g. We 
therefore postulate the following scaling hypothesis 



V? nt (V,L) = f(L/aV)) 



(13) 



where £(V) is the correlation length of the MIT. In the 
insulating phase (V > V c ) the conductance of a sample 
with length L behaves as g(L) ~ exp(— L/£) due to the 
exponential localization of the eigenstates, and therefore 
we have g{L{) < g{L 2 ) for L\ > L 2 . Based on Eq. ( p^ ) 
we expect the same behavior for V® nt i.e. for every finite 
L 1 > L 2 we must have 7^(7, Li) < V° nt (V,L 2 ). On 
the other hand, in the metallic regime (V < V c ) we have 
that g(L) = DL and therefore we expect from Eq. (|T^ ) 
V? at (V,Li) > P? nt (V,L 2 ). Thus, the critical point is the 
one at which the size effect changes its sign, or in other 
words, the point where all curves P® nt (V,L) for various 
L cross. These considerations are confirmed in Fig. ||a, 
where we report as an example the evolution of V^ nt (V) 
for different L using the box distribution. From our data 
we determine the critical point to be V = V c = 16.5 ±0.5 
in agreement with other calculations [p|Jl0[. In order to 
verify the scaling hypothesis (|l3| ) we report in Fig. ^Jd 
the same data as a function of the scaling ratio L/£. All 
points collapse on two separate branches for V < V c and 
V > V c . We consider this to be a good confirmation of 
the scaling hypothesis. 

We now turn to the analysis of the delay time statis- 
tics V(f). Equation (^]) provides us with a simple under- 
standing of the scaling f = tAM = rj < t > proposed 
in (|l|). In Fig. 4a we report the distribution of the log- 
arithm of the rescaled delay times ■p(ln(f)) p5[ for the 
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three different distributions V(V n ) of the random poten- 
tial and for various sample sizes L. In all cases the distri- 
bution functions lie one on top of each other, indicating 
that their shape is universal. Together with our previous 
result for V(T) this indicates that at the MIT only one 
fundamental energy scale, the mean level spacing A, is 
required to characterize the statistical properties of the 
5-matrix. 
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FIG. 3. (a) Universal behavior of V(f) at the MIT (re- 
ported here as P(ln(f))) for various sample sizes L and po- 
tential distributions, (b) The integrated distribution "Pint ("?) 
asymptotically decays as "Pint (f ) ~ f 15 (dashed line). 

We now discuss the general features of V(f). For f Cl 
the behavior of V(f ) can be deduced from Fig. |^a, which 
clearly shows that ln'P(lnf) decreases faster than Inf. 
This guarantees that V(f —> 0) — > since V(f — > 0) = 
c ^ implies ln"P(lnf) = ln(c) + hi(f). Thus, a gap is 
formed at the origin of V(f). We point out that a similar 
behavior is predicted for the conductance distribution at 
the MIT by the e-expansion in field theory pl| a nd has 
been verified by recent numerical calculations [ 10| . 

In Fig. [|b we present V- m t(r) in a double logarithmic 
plot. Again all curves coincide. Moreover the tails show a 
power-law decay 7 3 i n t('?) ~ with 7 w 1.5 ±0.05 given 
by a least square fit. The fact that 7 > 2 indicates that a 
mean delay time can be defined at the MIT in contrast to 
the localized regime where we have divergence of the first 
moment of V(t) ||. This is yet another manifestation of 
the fact that at MIT the conductance has a well defined 
value g c . We note the divergence of higher moments of 
V(t ) and point out that a similar behavior is predicted 
for the higher cumulants of the conductance distribution 
at the MIT from field theoretical calculations |ll| . 
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